Creating perturbations from a decaying field during infiation 
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Typically the fluctuations generated from a decaying field during inflation do not contribute to 
the large scale structures. In this paper we provide an example where it is possible for a field which 
slow rolls and then decays during inflation to create all the matter perturbations with a slightly 
red-tilted spectral index, with no isocurvature perturbations, and with a possibility of a departure 
from Gaussian fluctuations. 
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Primordial inflation is one of the most successful 
paradigms for the standard model cosmology which has 
many observational consequences pj . In spite of this suc- 
cess, understanding the origin of inflation is still very 
challenging [2 . In the simplest version it is assumed 
that a light scalar field is slowly rolling down the po- 
tential which produces almost scale invariant Gaussian 
fluctuations [3]. 

However there could be departures from this simplest 
paradigm. In principle there could also be many scalar 
fields which might participate during inflation [3], out of 
which some could be heavy and they might decouple from 
the dynamics early on, while leaving the lightest field to 
slow roll to end inflation, or some could be ultra light 
which might even decay after the end of inflation, such 
as in the case of a curvaton scenario [Sj . 

One such physical possibility which we wish to explore 
here is the fate of a heavy scalar field during inflation. 
Since the field is very heavy compared to the Hubble 
expansion rate, one might expect it to settle down in 
its potential within one Hubble time. Its perturbations 
would be damped and as a result it would not leave any 
observable imprint on the large scales, see for instance [B]. 

However this might not be the case if the heavy field 
lies initially on a plateau that is flat enough to accommo- 
date slow roll, such as in the vicinity of a saddle or inflec- 
tion point, very similar to the cases of inflation discussed 
in Refs. j7j. If such a field rolls off the plateau during the 
last 50 — 60 e-folds, then it is possible to seed the pri- 
mordial fluctuations for the observable scales which will 
be dominated by the fluctuations of the heavy field while 
making the inflaton's own fluctuations redundant. 

In order to probe such a scenario, we need to under- 
stand the dynamics of both the fields, i.e. the heavy and 
the inflaton. Let us consider an inflaton (p with poten- 
tial V{(j)), and another field a with potential t/(cr) added 
upon it, so the total potential is 



Vtot{^,a) = Vi(f>) + Uia). 



(1) 



We consider the case where the inflaton (j) leads inflation, 
so V > U always holds, whereas perturbations of 4> are 
negligible compared to that of a. We also assume V{a) to 
have a plateau flat enough to accommodate slow roll. To 
simplify the calculations we may assume that the U{a) 
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FIG. 1. The schematic timelines of background and perturbed 
universes. The blue, red and green solid lines denote a, <j) and 
total energy densities of the background universe respectively. 
The dashed green lines show the total energy densities in the 
perturbed universes where the main contribution comes from 
fluctuations in the a fleld. After the Hubble exit, inflation is 
divided into two phases by the break of ct's slow roll condition. 



potential turns down very sharply at the plateau edge. 

Because of the plateau in U{a), we can split inflation 
into two distinct phases with three critical points, whose 
timelines are shown in FIG. [T] During the flrst phase, 
both the fields are slowly rolling as a remains on the 
plateau. This is terminated by the break of the second 
order slow roll condition of cr field as it reaches the edge 
of the plateau. Then a rolls quickly to its minimum of 
the potential and oscillates, which we regard here as the 
second phase. During this phase a oscillates perhaps 
many times within one Hubble time and its quanta are 
quickly redshifted away by inflation. Several e-folds af- 
ter cr ends slow roll, it is already diluted enough so the 
universe enters a single fleld inflation. 

As shown in FIG. [T] perturbations in cr fleld is actu- 
ally responsible for a fluctuating boundary between the 
phases, generating different Hubble rates around the time 
a ends slow roll. This imparts fluctuations in the number 
of e-folds: 



SN = N^Sa^, 



(2) 



where — dN/da. The power spectrum for curvature 
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perturbation can be given by: 



Pc = Ps 



N 



(3) 



in which Psa, ~ HJ^/^tt'^ for any sufficiently flat poten- 
tial during inflation, with iJ* being the Hubble rate of 
expansion when the relevant modes for observations are 
leaving the Hubble patch. 

Exact analytical calculations exist for two-field infla- 
tions [S], but we can estimate Pq easily by using the 
evolution of the inflaton as a clock. As shown by the 
dashed green lines in FIG.jlj we see that the fluctuations 
in a lead to different locations for the end of slow roll 
condition in perturbed/background universes. As an es- 
timation, we may assume that once a is off the plateau 
it doesn't contribute at all to the energy density. This 
simplifies pa to almost a step function instead of the one 
shown in FIG.jl] Then the initial perturbation Jcr* moves 
the location of "c", (which indicates when a ends slow 
roll) by the amount of time 
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or by the amount of inflaton clock time: 

V' 



= SrSt 



(4) 



(5) 



From eq. ([9]), we see that the result we have estimated 
for this scenario is very simple. The spectrum is solely 
governed by the dynamics of U{<j) instead of V {(/)), and 
the spectral tilt of the curvature perturbations is domi- 
nated by: 



Us - I 



dlnR 



C 



Hdt 



-2e* 



where 



d_ J_ 
dtWl 



(10) 



(11) 



is the total flrst order slow roll parameter. From eq. ( 10 ) 



we find the the red-tilt is favored by current observations, 
i.e. the central value of the spectral tilt is given by: Ug w 
0.96 [T]. Furthermore, since inflation goes on, — 1 will 
be further negative, the running in the spectral tilt will 
also be slightly negative, which is also favored by the 
current observations. 

After providing the readers with this heuristic argu- 
ment, let us now proceed with a detailed calculation 
which will also yield the the bispectrum of curvature 
perturbations. Starting from the flrst phase, the slow 
roll approximations give the equations of motion for both 
flelds 



During this time in a perturbed universe a would still 
contribute to the total energy density, while it would not 
in the background universe. This brings about the dif- 
ference in the Hubble rate by: 



tt2 

pert 



u 



v + u 



(6) 



From the slow roll approximations we get the numbers 
of e-folds of inflation that are generated during this infla- 
ton clock time S(j)c, for the background and the perturbed 
universes 



N 



3fP 

'V7 



pert 



"^^pert 
V' 



(7) 



Since the only difference between the perturbed and the 
background universes' evolutions comes from the period 
54>cj we arrive at[^ 



5N EE Npert " N 



so that 



(8) 



(9) 



^ We have approximated Uc ~ Ut for a flat enough plateau. 



d(/)/dt = ^V'/ZH, Aa/dt = -U' /3H , (12) 

whose ratio after integration is 



V 



da 

. U' 



(13) 



When being perturbed by the initial Sa^,, because ct's 
plateau edge is sharp it will always end slow roll at almost 
the same position, regardless of the perturbation. This 
allows us to neglect Sac to get the perturbed version of 



eq. (131 



Ve 



Sa, 



= 0. 



(14) 



Then it is easy to write down the perturbed number of 
e-folds during the first phase, for a fiat enough plateau 



SNi = - 



8Tr{Vc + Uc) 



SiriVc + Uc) 



M^Vc 



M^U, 



Sa^, 



(15) 



After a ends slow roll, it may either oscillate around its 
vacuum with a frequency much higher than the Hubble 
rate, or it may decay instantly. Irrespective of these two 
scenarios we will have more or less a constant overall 
equation of state, w, for the a oscillations or its decay 
products. In either case the decay process is generally 
unimportant in this scenario as everything will be diluted 
away by the follow-up infiation. 
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Let us consider the second phase when it starts as a 
ends slow roll. From the Friedmann equation and the 
equation of motion for the slow roll field we get 

where 0' = d0/dA^ = —d(f>/Hdt. This equation is a first 
order differential equation between (f> and N, which has 
an exact solution for the number of e-folds of the second 
phase: 



N2 = n{(f)c,(t>e) + 



f 



3(f + «;) 



In f 



r(l — a) 

I - r 



where 



8nV 



(17) 



(18) 



is the number of e-folds for the part of a single field infla- 
tion with (f) varying from 0i to 02 , as if the energy density 



from (T did not exist. The second term in eq. (17) is the 



additional number of e-folds provided by the presence of 
the oscillating a or its decay products, however negligible 
the contribution is. Also 



(19) 



and a <C 1 has the order of the inflaton's first order slow 
roll parameteiF] 



After calculation, the perturbation of N2 from eq. (17) 
can be expressed as: 



6N2 = - 



1 - r 8n{Vc + Uc) 



1 



5a^ 



(21) 



which only differs from eq. ( 15 ) slightly. The dependence 



on the equation of state w is fully encoded in a. The total 
difference in the number of e-folds of inflation comes from 



combining eq. ( 15 ) and eq. ( 21 ) , during which most terms 



cancel because of the fluctuating boundary "c" , yielding: 



SN = SNi + SN2 = NJa^ 



where 



37r(l - a)Uc 



(1 - ar)MlU^' M^U^' 



(22) 



(23) 



Since a is much smaller than unity and « U we 
confirm our previous estimation result in eq. ([8|. Also, 



^ If the inflaton <f> remains almost constant we will get q — > 0. The 
exact definition of a is given by: 

2A{1 + w)ttVc e-3(i+™)"(<#'c.« 



a = 1 



Ml 



■f 



V 



(20) 



we get the power spectrum of curvature perturbation (ac- 
cording to eq. (I3I) 



16(1 - a)2i7,2[/^2 iQH.,^Uc 



(1 



ar; 



(24) 



From the exact expression of eq. (24), we also get the 



exact spectral tilt for the curvature perturbation as: 



l.^ = -2..+2,. 



(25) 



where ?7^cr* = (Mp/87r)[/*"/(y* + C/*) is the second order 
slow roll parameter for a. 



With the exact expression eq. ( 23 ) , we can now calcu- 
late the bispectrum of the curvature perturbations. By 
taking the derivative d/da^, on both sides of the exact eq. 
(23), we can define and calculate N^ra = dNg-fda^. After 
some work, we get the strength of the local bispectrum 
for the curvature perturbation as: 



/i 



NL 



6 



6r ' 



(26) 



where /3 is mostly a function of a, r, w and e^c = 
MpVc'^ /WttVc^ the slow roll parameter for the inflaton 
itself when a ends slow roll. 

Let us now consider a very special case when ?7(cr) is 
very subdominant so V{4>) 3> U{cr), and a potential is 



flat enough on the plateau. We can then simplify eq. ( 24 ) 
to 



,1/2 



4 



and for eq. (26), /? reduces to 



3(1 



1 - a 



2e0c < 1 . 



(27) 



(28) 



Since the numerator and the denominator are both much 
smaller than 1 in eq. (26), whether we get a signifi- 



cant non-Gaussianity then depends on which one is the 
smaller. 

Unlike the power spectrum of curvature perturbations 
which is determined by U{(j), the value of (3 is decided 
solely by the inflaton potential V{(j)) and the equation of 
state w. We can predict the bispectrum on a per inflation 
model basis by calculating /?. As an illustration, we take 
V{(j)) = rn^cj)^, and the assumption that a instantly de- 
cays into radiation once it ends slow roll. As displayed in 
FIG.[2] we obtain /3 < while \ f3\ ~ 0(10-*) - 0(10-2), 
depending on the e-folding when a ends slow roll. 

Furthermore, since the decay products of a dilute away 
during inflation, there will be no relic isocurvature per- 
turbations. In this respect our paradigm is very different 
from the curvaton scenario as discussed in Refs. One 
of the challenges for a curvaton paradigm is to ensure 
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FIG. 2. The /3 parameter for V{4>) = m?4>'^ model when a is 
subdominant. The value of /3 depends on A*', the remaining 
number of e-folds of inflation after a ends slow roll. Note here 
= indicates end of inflation and A'^ ~ 50 is the pivot scale. 



that both curvaton and inflaton decay products thermal- 
ize which is highly non-trivial. If the curvaton does not 
dominate the energy density, one has to ensure this hap- 
pens otherwise large isocurvature perturbations can be 
produced [9j . But in this current scenario the origin and 
the decay products of a decaying field do not lead to any 
imprint in the thermal history of the universe other than 
the curvature perturbations. 

On another account our model predictions differ from 
the curvaton scenario [S]. The curvaton mechanism nat- 
urally predicts almost scale invariant perturbations, but 
in our case we naturally obtain red-tilted spectral index. 



see eq. (25 1, in accord with the current observations [T]. 

The flatness of a field can be ensured within particle 
and string theory [2] . However now the onus is to embed 
V{(t>) within the observable sector such that the infla- 
ton decay products directly excite the Standard Model 
quarks and leptons. Once this is ensured it is sufficient 
to create all observed matter and perturbations. In this 
respect, the decaying field might even originate from the 
ubiquitous hidden sectors available in beyond the Stan- 
dard Model physics. 

Before we conclude our results, let us briefly mention 
that decaying scalar field during inflation can produce 
significantly large /nl depending on the value of r. For 
an example if r sa 10^'* for the case when V{(t>) — m 



with |/3| ~ 0(10^), the non-Gaussianity parameter could 
be of order /nl ^ C(10). If some mode leaves the Hub- 
ble patch right before a ends slow roll, the Jml could be 
made even larger since the suppression factor P ~ 
However the running of the spectral tilt could also be 
significant enough to be ruled out by the current obser- 
vations. Some of these issues can be investigated in future 
publications. 

Let us conclude by mentioning that we have provided 
a simple paradigm where a decaying scalar field during 
inflation can indeed seed red-tilted spectrum with the 
desired amplitude of perturbations which can give rise 
to small or large f^L depending on the ratio of r and 
the potential of inflation at the time when the decaying 
field leaves the slow-roll approximation. 



AM is supported by the Lancaster-Manchester- 
Sheffield Consortium for Fundamental Physics under 
STFC grant ST/J000418/1. 



[1] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. 

Suppl. 192, 18 (2011) arXiv:1001.4538 [astro-ph.CO]]. 
[2] A. Mazumdar and J. Rocher, Phys. Rept. 497, 85 (2011) 

arXiv:1001.0993 [hep-ph]]. 
[3] V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, 

Phys. Rept. 215, 203 (1992). 
[4] A. R. Liddle, A. Mazu mdar and F. E. Sch unck, Phys. Rev. 

D 58, 061301 (1998) astro-ph/9804177'. 
[5] D. H. Lyth and D. Wands, Phys. Lett. B 524, 5 (2002) 
hep-ph/0110002 . D. H. Lyth, C. Ungarelh and D. Wands, 

Phys. Rev. D 67, 023503 (2003) astro-ph/0208055 . 

K. Enqvist and M. S. Sloth, Nucl. Phys. B 626, 395 (2002) 
hep-ph/0109214 . T. Moroi and T. Takahashi, Phys. Lett. 

B 522, 215 (2001) [Erratum-ibid. B 539, 303 (2002)] hep^ 

ph/0110096. 

[6] A. RLiddle and A. Mazumdar, Phys. Rev. D 61, 123507 

(2000) astro-ph/9912349 . 
[7] R. AUahverdi, K. Enqvist, J. Garcia-Bellido and 
A. M azumdar, Phys. Rev. Lett. 97, 191304 (2006) [hep^ 
I ph/0605035 . R. AUahverdi, K. Enqvist, J. Garcia-Bellido, 



A. Jokinen and A. Mazumdar, JCAP 0706, 019 (2007) 
hep-ph/0610134. 
[8] A. Mazumdar and L. -F. Wang, JCAP 1209, 005 (2012) 

arXiv: 1203.3558 [astro-ph.CO]]. 
[9] A. Mazumdar an d S. Nadathu r, Phys. Rev. Lett. 108, 
111302 (2012) arXiv: 1107.40781 [hep-ph]]. 



